Quantum phase transition in capacitively coupled double quantum dots 
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We investigate two equivalent, capacitively coupled semiconducting quantum dots, each coupled 
to its own lead, in a regime where there are two electrons on the double dot. With increasing interdot 
coupling a rich range of behavior is uncovered: first a crossover from spin- to charge-Kondo physics, 
via an intermediate 57/(4) state with entangled spin and charge degrees of freedom; followed by a 
quantum phase transition of Kosterlitz-Thouless type to a non-Fermi liquid 'charge-ordered' phase 
with finite residual entropy and anomalous transport properties. Physical arguments and numerical 
renormalization group methods are employed to obtain a detailed understanding of the problem. 
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Introduction. - Semiconducting quantum dots provide 
[1] a beautifully direct, tunable mesoscopic realization 
of a classic paradigm in many-body theory: the spin- 
Kondo effect [2] , wherein a single spin in an odd-electron 
dot is quenched by coupling to the conduction electrons 
of a metallic lead. Recent advances in nanofabrication 
techniques now also permit the controlled construction 
of coupled quantum dot systems, the simplest being dou- 
ble dot (DD) devices. Central to the design of circuits 
for logic and quantum information processing, and widely 
studied both theoretically [3-9] and experimentally [10- 
15], spin and orbital degrees of freedom are now relevant, 
leading to the possibility of creating novel correlated elec- 
tron states. Recently for example a symmetrical, capac- 
itively coupled semiconducting DD has been studied [5] 
in a regime with a single electron (n = 1) on the DD, and 
the lowest energy states (nt, tir) = (1, 0) and (0, 1) near 
degenerate. The low-energy physics, which determines 
the conductance at small bias, was shown [5] to be gov- 
erned by a fixed point with SU (4) symmetry, leading to 
an unusual strongly correlated Fermi liquid state where 
the spin and orbital degrees of freedom are entangled. 

In this paper we study a capacitively coupled, sym- 
metrical semiconducting DD system, but now in a regime 
with two electrons on the DD such that (1, 1), (2, 0) and 
(0, 2) are the relevant low energy states. As shown be- 
low, the associated physics is both rich and qualitatively 
distinct from the n = 1 sector: on increasing the ra- 
tio U'/U of inter- and intra-dot coupling strengths, we 
find that the system first evolves continuously from an 
SU(2) x SU(2) spin-Kondo state where the dot spins are 
in essence separately quenched, to an SU(4) Kondo state 
with entangled charge and spin degrees of freedom when 
U'/U = 1. Thereafter, for a tiny increase in U'/U, there 
is then a smooth crossover to a novel charge-Kondo state; 
followed by suppression of charge-pseudospin tunneling, 
manifest in collapse of the associated Kondo scale and a 
Kosterlitz-Thouless (KT) quantum phase transition to a 
doubly degenerate charge-ordered state — a non-Fermi 
liquid phase with In 2 entropy and anomalous low-energy 



transport and thermodynamic properties. Detailed re- 
sults for this diverse range of behavior are obtained us- 
ing the numerical renormalization group (NRG) method 
[16, 17], preceded by simple physical arguments that en- 
able the essential physics to be understood. 

Model and physical picture. - We consider two equiv- 
alent, capacitively coupled semiconducting (single- level) 
dots, each coupled to its own lead. The Anderson- type 
Hamiltonian is H = Hq + H\z + Hr>, where Hq = 
Yli k a £ k a ticr a ki<T re f ers to the leads (i = L/R) and 
Hy = J^i k a ^( a L<T C ier + h.c) to the lead-dot couplings. 
Hr> describes the isolated dots, 



H d = }^ (ehj + Uh^hji) + U'h L h R 



(1) 



i=L,R 



with hi = ^Zo^io = So- c \a c ia- U denotes the intradot 
Coulomb interaction, U' the interdot (capacitive) cou- 
pling. In the isolated DD, increasing |e| = — e (via suit- 
able gate voltages) generates the usual Coulomb block- 
ade staircase. For < |e| < min({7, U') the ground 
state occupancy is n = 1, with degenerate configura- 
tions (riL,nR) — (l,0)/(0, 1) and a ln4 residual entropy 
(fog = 1) [3-5]. Coupling to the leads quenches this en- 
tropy, and the strongly correlated effective low-energy 
model is SU(A) Kondo [4, 5]. The underlying physics 
here is rich, including spin-filtering arising from the con- 
tinuous crossover to the SU (2) orbital Kondo effect in a 
strong magnetic field [5]. But no quantum phase transi- 
tion occurs in this n = 1 sector. 

We consider by contrast the n — 2 domain of the 
Coulomb staircase, arising for mm(U,U') < |e| < U'+ 
m&x(U,U'). Two sets of configurations then domi- 
nate, according to whether U' ^ ^ : the 4-fold spin- 
degenerate states (nL,nu) — (1,1), and the degenerate 
pair (2,0)/(0,2), with DD energy difference £(2,0) - 
E(l, 1) = U- U'. The DD ground state is thus (1, 1) for 
U' < U, and (2, 0)/(0, 2) for U 1 > U; all six states are de- 
generate at U' — U where the model has SU (4) symme- 
try. We first give physical arguments for the evolution of 
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the coupled DD-lead system with increasing U'\ focusing 
on the strongly correlated regime of U/T 3> 1 (r = nV 2 p 
with p the lead DoS). Here an effective low-energy Hamil- 
tonian may be obtained from second order perturbation 
theory (PT) in the lead-dot tunneling V = Vl = Vr 
(with V L m denoting coupling to the L/R lead) [18]. 

For U' = the dots are fully decoupled. Only (1, 1) 
states are relevant. The effective model is obviously two 
uncoupled spin-^ Kondo models. The spin entropy is 

quenched at the normal Kondo scale T K U ^ 2 \ leading to a 
local singlet ground state ( l SU(2) x SU(2) spin-Kondo'). 
For U' > U by contrast the (2,0)/(0,2) DD states dom- 
inate, and as shown below the ground state is a doubly 
degenerate 'charge ordered' (CO) state with In 2 entropy. 
Continuity then implies a quantum phase transition at 
some critical U' c . As discussed below, for U' = U the 
effective model is SU(4) Kondo (in the n = 2 sector) 
[18], with entangled spin/charge degrees of freedom but 
a singlet ground state with a larger Kondo scale T K U ^; 
and is connected continuously to the SU (2) x SU(2) spin- 
Kondo state arising as U' — > 0. We thus expect U' c > U. 

Hence consider increasing U' above U. Since the con- 
figurations (1, 1), (2, 0), (0, 2) are degenerate for U' = U 
this full n = 2 manifold must be retained for U' ~ U. 
Virtual excitations to excited states are eliminated via 
PT, and divide in two classes [18]: (a) V% or Vj^ processes, 
involving tunneling to one lead alone. Any configuration 
connects to itself via such, eg (1,1) (1,1) under V/j via 
excited states (1,0) or (1,2). (b) VlVr processes. These 
necessarily connect different manifold configurations; the 
full set is clearly (2, 0) (1, 1) and (0, 2) <-> (1, 1), there 
being no direct coupling between (2, 0) and (0, 2). As U' 
increases above U charge and spin states begin to sepa- 
rate: the degenerate charge pair (2, 0)/(0, 2), components 
of an effective charge pseudospin, lie lower in energy by 
V — U than the (1,1) spin states. For sufficiently small 
U' - U > 0, tunneling between the (2,0)/(0,2) states 
can however still arise, and quench the charge pseudospin 
(and hence entropy) , producing thereby a non-degenerate 
charge-Kondo state. But, as above, this tunneling is 
not direct, being mediated by the higher energy (1,1) 
states. We thus expect the associated Kondo scale to be 
diminished compared to T K U ^ (the 'stabilization' due 
to dot-lead coupling at the SU(A) point) and to decrease 
as U' increases; and moreover that the quenching will 
cease to be viable when the relative energy U' — U of the 
(1,1) states exceeds roughly T K U ^ 4 \ leading to a quan- 
tum phase transition to the degenerate CO phase when 
U' c — U~ T K U< " 4 \ Since the latter is exponentially small 
for strong correlations, this implies a critical U' c exponen- 
tially close to U (as confirmed by NRG below). 

These arguments extend readily to U' < U, but now 
the circumstances differ. Spin/charge degrees of freedom 
do separate on decreasing U' from U, (1,1) states now 
lying lower by U - U' than (2,0)/(0,2). But since the 
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FIG. 1: (a) Phase diagram in (U' , (7)-plane; the SU(4) line 
U' = U is also shown (dotted), (b) For U = 7 in the SC phase, 
ln(TW/r) vs [U' c — U']^ 1 ^ 2 close to the transition, showing the 
exponential vanishing of Tk- Inset: ln(TW/r) vs U' . 



(1, 1) spin-states connect directly to themselves under Vj^ 
or V£ (as above), quenching of the spin entropy is not 
inhibited and no transition occurs. Instead a continuous 
crossover from SU(4) Kondo to separable SU (2) x SU{2) 



spin-Kondo is expected, for U — U' w T 
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Results. - The physical picture is thus clear, and the 
transition to the degenerate CO phase occurs in the vicin- 
ity of the SU(4) point U 1 = U [18]. We now present 
NRG results for the DD Anderson model, choosing the 
midpoint of the n = 2 domain, |e| = U/2 + U'. This case 
is particle-hole (p-h) symmetric, but representative. The 
physics is wholly robust to departure from p-h symmetry. 

The low temperature behavior of the model is gov- 
erned by two classes of stable fixed points (FP), corre- 
sponding to the two zero temperature phases. The first, 
a strong coupling (SC) fixed point, describes all the sin- 
glet ground states; and is reached (as T — > or NRG 
iteration number N — ► oo) for all U' < U' c . The cor- 
responding FP Hamiltonian is simply a doubled version 
(SU(2) x SU{2)) of that well known for the spin- 1/2 An- 
derson model [17]. The dot spins and hence entropy are 
thus quenched at T = 0, the system being a Fermi liq- 
uid and characterized by a Kondo scale denoted generi- 
cally as Tk- The second, reached for all U' > U' c , is a 
line (i.e. a one parameter family) of charge ordered (CO) 
FP. The generic FP Hamiltonian corresponds to setting 
r = and U' = oo. The DD and leads are then decou- 
pled, but the FP has internal structure reflecting broken 
symmetry, since dot-states occur in the degenerate pair 
(nL,nji) = (2,0)/(0,2) (whence ln2 residual entropy). 
The line of FP is obtained by supplementing the free lead 
Hamiltonians by potential scattering correlated to dot oc- 
cupancy, of form H K = Jf£ ij<T £ k>k , aL^k'ia^ -1 ) ( c f 
[18]). The actual value of K is obtained numerically by 
matching to NRG energy levels. 

A comparison of NRG energy level flows for large iter- 
ation number, with the characteristic energy level struc- 
ture for the two FP, enables the phase diagram to be 
found; as shown in Fig. 1(a) in the (U' = U'/ttT,U = 
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FIG. 2: S(T) vs T/r for U = 7 ({% ~ 7.046). (a) C/' = 6 (SC, 
'uncoupled spin-Kondo') and 8 (CO); (b) £/' = 6.9 (crossover 
to SI/(4)); (c) t/' = 7 (517(4)); (d) 17' = 7.03 (charge-Kondo). 
Dotted lines show In 16, In 6, In 4 and In 2, associated with the 
FO, LM sc/(4) , LM sc/(2) and CO fixed points respectively. 



U/ttF) -plane [19]. The transition is seen to occur for all 
U > on increasing the interdot U' . Consistent with 
the physical arguments above, for U 3> 1 the critical U' c 
indeed lies exponentially close to the SU(4) line V = 
U (specifically we find (U'JU - 1) ~ 2T K U(A) /TV 1 ' 2 ). 
Fig. 1(b) (inset) shows the [/'-evolution of Tk [19] in the 
SC phase, for a typical strongly correlated U = 7. It 
is seen to depart little from its U' = value T K U ^ oc 
rf7 1 /2 ex p(_i/ /9 j) ( w i t h P J = 8/n 2 U) [17], indicative of 
spin-Kondo physics, until very close to V = U where 
it increases rapidly to X K U ^ oc rt/ 3 / 4 exp(-l/2pJ); le 
while exponentially small, y^^ 4 ) oc [j 1 ^^ 2 )] 1 / 2 shows a 
strong relative enhancement at the 5E/(4) point [4, 5]. 
However on increasing U' above U and entering the 
charge-Kondo regime, Tk is seen to drop rapidly and 
vanishes as the SC— >CO transition is approached. The 
transition is of KT type, consistent with the line of CO 
FP for U' > U' c , and (from NRG energy level flows) no 
evidence for a separate critical FP; further evidenced by 
[20] the U' U' c - behavior T K oc exp(-a/[E/£ - l>] 1/2 j, 
demonstrated in Fig. 1(b) for U = 7. This behavior is 
generic, even for U = 0; here the non-interacting SU (4) 
{U 1 = 0) scale T K U{i) ~ T, and U' c -U ss T K U{4) implies 
U' c ~ 0(1), as indeed found (Fig. 1(a)). 

In addition to the two stable (low-temperature) FP, 
three unstable FP play an important role at finite-T, and 
are seen clearly in the impurity entropy 5 (= 5i mp ): (i) 
free orbital (FO) [17], corresponding to T = = U = U', 
with In 16 entropy. This is just the high-T limit of S(T), 
reached in all cases for non-universal T w max(?7, U'). 

(ii) SU(4) local moment (LM SU ^), corresponding to 
r = and U = oo = U' , with associated entropy In 6; and 

(iii) SU(2) x SU(2) local moment (LM SU ^), V = = U' 
and U = oo, with ln2 2 entropy. Fig. 2 shows S(T) vs 
T/r for U = 7 (l> ~ 7.046). Fig. 2(a) illustrates the 
behavior 'deep' in the CO and SC phases. In the former 
S(T) simply crosses directly from its In 2 (CO) residual 
value to In 16 (FO) on the scale T ~ U'\ while in the 
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FIG. 3: Transmission nTDi(u}) vs co/T, for [7 = 7 and (top 
to bottom) U' = 7, 7.03, 7.044 (SC) and 7.048, 7.1 (CO). 



latter, consistent with the effective underlying spin-1/2 
Kondo physics, there is first a crossover from 5(0) = 
(SC) to In 4 (LM SU W) for T ~ T K (^ T K U{2) ). Figure 
2(b), for U' = 6.9 < U, illustrates the crossover from 
effective uncoupled SU(2) Kondo to SU{4). Here 5(T) 
increases in a two-stage fashion, first to In 4 (LM SU ^) 
for T - T K and then In 6 (LM SU ^) for T ~ [U - U'] 
(= E(2, 0) — E(l, 1)). This behavior is naturally absent at 
the SU(4) point, Fig. 2(c): 5(T) crosses directly to In 6 
for T ~ T K U< " 4 \ In the charge-Kondo regime Fig. 2(d) 
(for U' = 7.03) 5(T) again shows the two-stage behav- 
ior typical of a KT transition [9, 21], but here it first 
crosses from (SC) to In 2 (CO) for T ~ Tk and then 
to In 6 (LM SU ^) for T ~ [U' - U] = £(1,1) - £(2,0); 
consistent with the physical discussion given above. 

The physics discussed above naturally shows up also 
in various thermodynamic susceptibilities. For example, 
as U' is increased past U into the charge-Kondo regime, 
we find the 'impurity' spin susceptibility \s decreases 
monotonically; but its charge pseudospin analogue the 
staggered charge susceptibility \~ , given by \c ~ 
diverges as U' — ► U' c reflecting the collapse of the charge- 
Kondo state and the quantum phase transition. For U' > 
U' c the (T = 0) Xc remains infinite, symptomatic of the 
broken symmetry CO phase, with Xc(T) oc 
0. Further details will be given in subsequent work. 

Finally and most importantly, the destruction of the 
Kondo effect as the SC^CO transition is approached is 
seen vividly in electronic transport, notably the transmis- 
sion coefficient T(lo) = ttTD^lo) with AM the T = 
local single-particle spectrum {Di(uo) = —ImGi{uj)/-n 
with Gi(t) = -i0{t)({c itT (t),c] a }) the retarded dot Green 
function). At the Fermi level in particular (w = 0), Tj(0) 
gives the linear differential conductance across one (ei- 
ther) dot in units of the conductance quantum 2e 2 /h 
[22]; and at finite, low bias voltage V, the equilibrium 
T{uo = eV) provides an approximation to the conduc- 
tance [22]. The low-energy behavior of Ti(uj) is shown 
in Fig. 3 for a range of U' spanning the transition. 
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For U' — U a relatively broad Kondo resonance char- 
acteristic of the SU{&) point is apparent, with width 
oc T K = T^ U ^ 4 \ On increasing U' into the charge-Kondo 
regime the Kondo resonance, now residing on top of an 
incoherent continuum, remains intact with T'j(O) = 1 
throughout the SC phase reflecting the unitarity limit 
[23]. But it narrows progressively as Tk diminishes, and 
asU' — > U' c — the Kondo resonance vanishes 'on the spot', 
such that for U' > U' c in the CO phase only the back- 
ground continuum remains. The linear conductance in 
particular thus drops abruptly at the transition. This ap- 
pears to be a general signature of an underlying KT tran- 
sition, it being found also for a multi-level small dot close 
to a singlet-triplet degeneracy point [21]; and in recent 
work [9] on two Ising-coupled Kondo impurities, onto 
which maps the problem of spinless, capacitively coupled 
metallic islands/large dots close to the degeneracy point 
between N and N + 1 electron states [8] . The non- Fermi 
liquid nature of the CO phase is also seen clearly here, be- 
cause Tj(0) = l/[l + (S I (0)/r)] in terms of the imaginary 
part of the dot self-energy; whence Tj(0) < 1 in the CO 
phase implies a non-zero T, l (ui = 0) and thus a non-FL 
state (as also manifest eg in anomalous exponents for the 
subleading T-dcpendcnce of thermodynamic properties, 
although the latter effects are quantitatively minor). 

We have considered an equivalent (L/R symmetric) 
DD system, with specific results shown for the p-h sym- 
metric case. For a DD device described by the effective 
model to be realizable, the physics should be suitably 
robust to breaking these symmetries. As mentioned pre- 
viously, that is entirely so for departure from p-h symme- 
try. L/R symmetry can be broken by detuning eg the dot 
levels, €r/l = t±St, or coupling to the leads, Tr ^ Tl- 
In that case it is readily shown that the SC FP remains 
stable, no new corrections to the FP being generated. 
For the CO FP by contrast, additional relevant pertur- 
bations arise. This FP is thus unstable, and flows in its 
vicinity ultimately cross over to the SC FP under renor- 
malization, with the crossover scale determined by St. 
We find however that eg for small but finite St <C T K U ^ , 
the thermodynamics above are essentially unaffected for 
T > St; and that the abrupt drop in the linear conduc- 
tance at the transition is simply replaced by a continuous 
but nonetheless sharp crossover over a U' interval on the 
order of T^ (4) . In that sense the physics described here 
is thus also robust to L/R symmetry breaking. 

Conclusion. - Motivated by extensive recent interest 
in capacitively coupled DD systems [3-15], we have an- 
alyzed a symmetrical, capacitively coupled semiconduct- 
ing DD in the two-electron regime, using the NRG tech- 
nique. We have shown that on increasing the intcrdot 
coupling the system evolves continuously through a pro- 
gression of Fermi liquid states from a purely spin-Kondo 
state, via the SU(4) point where charge and spin degrees 
of freedom arc wholly entangled, to a c/iarge-Kondo state 



with a quenched charge-pseudospin; before undergoing a 
quantum phase transition of the Kosterlitz-Thouless type 
to a non-Fermi liquid, doubly degenerate charge-ordered 
phase. This provides a striking example of the subtle and 
many-sided interplay between spin and charge degrees of 
freedom in small quantum dots. 
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